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Abstract 

Pairs of particles of definite total and relative angular momentum provide a 

natural description for a two dimensional electron gas in a strong magnetic 

field. Two body operators take a simple form when expressed in terms of 

pair creation and destruction operators. The pair formalism is applied to the 

study of edge waves excitations. For v = 1 the operators which create edge 

excitations are identified and the role the interaction potential plays in the 

long wavelength limit is clarified. This picture is claimed to describe also edge 

excitations on the u = 1/m Laughlin states. 
PACS. numbers 71.28, 71.10, 73.40 
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I. INTRODUCTION 



Contrary to the integer quantum Hall effect, which can be accounted for by a single 
particle description, the fractional effect arises as a result of condensation into a macroscopic 
collective ground state!. Much of the present understanding of the fractional quantum Hall 
effect (FQHE) is based on first quantized many electrons wave functions!'!. The strategy 
based on trial "variational" wave functions has the advantage of displaying in a very direct 
way the many body correlations between electrons. These correlations are induced by the 
pair potential acting on the electrons and this naturally leads to consider pair of electrons as 
the relevant degrees of freedom. Indeed the notion of particle pairs has been used by many 
authorslHzl in the study of the FQH effect. No systematic description of the system in terms of 
pairs and no detailed analysis of the the distribution of their angular momenta has been given 
up to now. The first part of this work is an attempt in this direction. In the symmetric gauge 
pair creation and annihilation operators are introduced in a second quantized formalism; 
their quantum numbers are the total angular momentum (TAM) and the relative angular 
momentum (RAM). All relevant two body operators in the lowest Landau level (LLL) are 
simply related to the distribution of pairs' TAM and RAM. In spite of the simple expression 
of the operators, the complex commutation relations among pairs, which deviate from perfect 
boson character, preclude a simple description of the Hilbert space. The first section is 
devoted to the derivation of the basic formulas, to the discussion of some simple examples 
and to a brief review of some well known results, easily recovered within the pair description. 

The second part deals with an application to a specific problem for which an interesting 
picture has been recently proposed on the basis of both analytical and numerical work!. 
This concerns the low lying excitations near incompressible states as the z/ = lorz/ = l/3 
states usually called edge waves. The pair picture provides a simple description of edge 
excitations. The main result is an operator relation that, under some conditions for the 
incompressible ground state, allows to identify the creation operators of edge excitations. 
The case v — 1, for which these conditions are satisfied, is treated in detail also relying on a 
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Hartree approximation. The same picture is suggested to hold also for the v = 1/m Laughlin 
states. The role the e~- e~ interaction plays in the edge wave dispersion relation is also 
clarified. In particular the results confirm the validity of the semiclassical approximation^ 
for edge states and suggest that the dispersion relation is asymptotically linear. Moreover it 
is shown that the contribution to the velocity of edge waves vanishes, in the thermodynamic 
limit, for any potential that decays faster than 1/rasr^ oo. 

II. PAIR DESCRIPTION OF THE FQHE IN THE SYMMETRIC GAUGE 

The Hilbert space of two dimensional charged particles in the x-y plane is split into 
Landau levels by a magnetic field B = B z z. In the extreme quantum limit (B z — > oo) 
all the particles are confined in the lowest level. The kinetic energy, reduced to the zero 
point motion, is an inessential constant so that the hamiltonian of the system contains only 
potential terms. Among these the dominant role is played by the particle-particle interaction 
potential. With this general motivation in mind we will concentrate in this section on the 
projection of two body operators on the LLL. 

Let us consider a general two body operator represented, in first quantization, by the 
function V{z\,z 2 ) (zj = Xj + iyj is the complex coordinate for particle j on the complex 



plane). Magnetic units, £ = ^hc/eH = 1, will be used throughout. The second quantized 
form of this operator in the LLL, in the symmetric gauge, is: 

i oo T 

^=oE Y.( T - u MV\T-s 1 s)4_ u c + u c s c T _ s (1) 

/ T=l u,s=0 

The operator c+ creates one electron in the LLL orbital (j) u {z) = (2tt2 u u\)~ 1 ^ 2 z u e~^ 2 ^ and 

c u is its hermitian conjugate. The conservation of TAM is explicit in equation (JTj) and the 

matrix element is given by: 

irp \ Iy(T,s,u) 

(T — u, u V T — s, s) — ; = (2) 

4vr 2 2 T J(T- s)\s\(T -u)\u\ 



with I v (T,s,u) = Jd 2 z 1 d 2 z 2 z 1 T - u z 2 u V(z ll z 2 )zl~ s ' z ^ e -(\^?+\^ 2 )/2 _ We consic i er i n w hat 
follows only operators V(zi,z 2 ) that are separable in the relative (£ = [z\ — z 2 )/2) and 



center of mass (Z = [z\ + z 2 )/2) coordinates: ViZ + £, Z — £) = w(Z) ■ v (£) In this case the 
integral in the matrix element can be performed0 as follows: 



T~S S T—U U / rp \ 

o=0 /3=0 7=0 o-=0 



a 



:-i) 



The integrals vanish whenever a + (3 ^ 7 + cr. Using a differential representation of the 
discrete delta function 5y = d^/^UM) and introducing the integer variable g = a + (3 = 
7 + a the sums on a, (3, 7 and a can be carried out and 

<?' 



%l + y) T - s (l-yy 



It(q)Iz(T,q) 



with Is(q) = 2 /rf 2 ^(0le| 29 e- |f|2 , I z {T,q) = J d 2 Z w(Z)\Z\ 2 ^e 



It is now possible to sum independently on u and s in equation ([!]), so that the second 
quantized form of the operator V becomes: 



00 T 

^EE« + (?)/r(?)' 

T=l g=0 



(3) 



The sum on u in equation ([!]), which involves the first integral and the pair of operators 
c t-u c u > defines the pair creation operator 



(4) 



u=0 



Similarly the sum on s yields the hermitian conjugate fx(q)- Requiring normalization of the 
state (<?)|0), some elementary algebra yields: 



b T {u,q) 
and 








G) 


2 T+1 



(l + z) T ~ u (l-z) u 



(5) 



2 = 



V T (q) 



Iz(T,q).I e (q) 



2Tt 2 {T-q)\q\ 

An orthogonality relation can be easily derived from equation 
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(6) 



T 

2J2 b T(u,q)b T (u,p) = 5 q , p (7) 

u=0 

since bT{u,q) = bx(q, u) an analogous relation holds for sums on q. This relation is of 
frequent use when q runs only on even (odd) values. In this case the sum can be extended 
to all q by inserting (1 ± (— 1) 9 )/2 in the sum. This yields 

4 b r(u, q)b T (s, q) = 5 UjS ± 5 UiT _ s (8) 
i 

the upper (lower) sign refers to even (odd) q. 

In the limit T — > oo the coefficients &t(w, q) are simply related to Hermite orthonormal 
functions <fi q {x) = ir~ (x) are Hermite polynomials): 

q 



b T U2/Tx + JT/2,q) = (2T)-4 



(9) 



Note that briT — u,q) = (— 1) 9 6t(w, q) so that, using the commutation relations of the 
operators c u and c+, we realize that for fermions fr{q) = for q even. Similarly, as a 
consequence of statistics, only even values of q occur if bose particles are concerned. This is 
confirmed by the symmetry of the first quantized wave function ijjT,q(Z, £) of the pair (T, q) 
under the exchange £ — > — £: 

ifr >q (Z, = V2 b T (% q)MZ + 0<t>T-u{Z - = 

7T—qcq 

- Z % p -\Z\V2-\tf/2 



2n^(T-q)\q\ 

This equation explicitly shows that the pair quantum numbers q and T represent the RAM 
and the TAM of the pair respectively. In a classical picture the pair's center of mass revolves 



round the center of the disk at a radius J (Z 2 ) = \JT — q + 1 while the two particles rotate 



around the center of mass Z on a circumference of radius J (£ 2 ) = y/q + 1. The average 



distance from the center of the disk of one of the two electrons is y (Z 2 + £ 2 ) = \/T + 2. 
This consideration is useful to identify the degrees of freedom relevant for the behaviour of 
the system in the bulk an at the boundary. A generalization to higher Landau levels of the 
expansion in pair wave functions in RAM and TAM was used by A.H.MacDonald et al.i. 



Let us now turn to the discussion of equation (H). In the case w(Z) = 1, Iz(T,q) = 
ir(T — q)\ and Vr(q) = e q is independent of T so that: 

£ '=H=f"(^)r^ (10) 

Note in particular that e q ~ v Cy/o) f° r q —> oo. The decomposition of pair operators in 
components of different RAM has already been introduced by Haldanei for the spherical 
geometry. In the disk geometry the same decomposition was studied by Trugman and 
Kivelson§ where a short range pair potential was expanded in powers of its range. An 
inversion formula for v(£) as a function of e q is obtained by multiplying equation (|10D by 
(— and summing over q: 



v(r) = 



n=0 



2£(-i; 



n 



L n (r 2 ). (11) 



q=0 

where L n (x) is the normalized Laguerre polynomial of order n. 

Examples of two body operators which can be expressed in the form (^j are: 

i) the Coulomb interaction V(z\, z 2 ) = l/\z\ — z<i\: 



V T ( q )=el = ^-^- 2 (12) 
y 1 q 2 (q\2i) 2 y ' 

note that for q ^> 1 e q oc l/^/q as in the classical case, 
ii) The hard core (HC) potential V(zi,z 2 ) = 27rV 2 5(z 1 — z 2 ): 

Vr(q) =e h q = 5 qA . 

iii) The pair correlation function g(r) corresponds to V(zi, z 2 ) = 5(r — Z\ + z 2 )\ 

VM = 9q (r) = ^- q e~ r2/4 . (13) 

Once the number of pairs with a definite RAM and TAM NT{q) = (fx (o , )/r( ( ?)) i s known 
on a given state for all T and q we are in a position to evaluate all relevant quantities. 
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The operators (<?) do not represent true bosonic particles. In fact the commutation 
rules for f^{q) has a residual term which contains a density excitation: 

= $T,RSq,p - (14) 

- 4 b r(u, q)b R (u, p)c^_ u c T _ u 

u 

The occurrence of such weird commutation relations makes the description of the Hilbert 
space in terms of pairs problematic. While only n/2 pairs are necessary to build a state, 
this contains much more pairs. Indeed it is easy to check that 

oo oo T i 

EE^(?) = EE c + u c u 4_ u c T _ u = -n(n - 1). (15) 
t=i q t=i«=o z 

This trivial sum rule simply states that n(n — l)/2 pairs can be made out of n particles. 

As an example of the pair's momentum distribution function it is easy to check that for 

the v = 1 state (all orbitals occupied up to n — 1) all pair states are occupied for T < n, i.e. 

N T (q) = 1 Vg, while 

n-l 

N T (q) = 2 b T(u,q) iorT>n (16) 

u=T-n+l 

Figure [I] shows the distribution of RAM N(q) = J2t ^t{q) for z/ = 1/3 and 6 electrons. Full 
dots refers to the true ground state of the Coulomb interaction while open squares to the 
state made of three pair creation operators with T = q = 15 acting on the vacuum. The 
latter state contains a lot of pairs with small RAM, including angular momentum one. The 
true ground state is instead characterized by a minimal number of pairs with the smallest 
value of q, since these give the largest contribution to the repulsive energy. In particular 
the Laughlin state (open dots) does not contain any pair with q = 1 so that it is the exact 
ground state for the HC interactionS'i. The dependence of N(q) on q has been analyzed 
numerically for system of n < 8 electrons. A fairly good scaling of the form N(q) = A q n — B q 
has been observed with A 3 = 0.7965, A 5 = 0.5856, Aj = 0.4659 and B q ~ 1 for all q. 

Equation ( |T3"D expresses a connection between the coefficients of different powers of r 
in g(r) and the number of pairs with a definite RAM. These coefficients have been studied 
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extensively by Yoshioka0, in rectangular geometry (Landau gauge), who has found that the 
coefficient of r 2 and of r 4 (which is non zero in this gauge) decreases by decreasing v and 
vanishes for v > 1/3. The same happens to the coefficient of r 6 and r 8 for v ~ 1/5. A similar 
result was derived by Trugman and Kivelsonil in the symmetric gauge. In the pair language 
the quantization in the FQHE occurs as a consequence of the successive elimination of all 
the pairs with the smallest RAM. 

III. EDGE WAVES IN THE QUANTUM HALL EFFECT 

We consider, in this section, the spectrum of low lying excitations on incompressible 
quantum Hall states in disk geometry. These excitations have been called edge waves (EW) 
because they involve density fluctuations of the two dimensional electron gas at the boundary 
of the system. A review on the subject can be found in a recent paper of X.G.Weni where 
a general theory for edge excitations is discussed. The starting point of Wen's theory is a 
classical hydrodynamical approach where coordinates and canonical momenta yield, upon 
quantization, the creation operators of the edge modes. The outcome of this approach is 
a free phonon theory. This is obvious if only the one body confining potential, coming 
from electron-background interaction, is considered. This picture however holds even in the 
presence of the e~- e~ interaction. A strong evidence of this has been given by M.Stone 
et al.i which have analyzed the energy spectrum due to the pair interaction using exact 
diagonalization for systems of up to 400 particles near v — 1. The ground state \ipo(n)) 
for v = 1 and n electrons is in first quantization (apart from the gaussian factors) the 
Vandermonde determinant ^ (zi, . . . , z n ) = Yii<j( z i — z j) whose total angular momentum is 
L Q = ^n(n — 1). In the sector of total angular momentum L = L + M the energy spectrum 
reduces, with excellent accuracy!, to 



E = E -Y j nA (17) 

k 

where rik are (integer) bosonic occupation numbers such that J2k n kk = M and Qk > 
are single particle energies. Note that, since only the interaction potential is considered, 
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the contribution to the EW spectrum is negative as a consequence of a loss of repulsive 
energy. If strictly only the LLL orbitals are considered, the Hilbert space of the system with 
L = L Q + M is spanned, in first quantization, by the wave functions obtained by multiplying 
\l/ by symmetric polynomials of degree M. These in turn can be expressed in a unique 
way0 as polynomials of power sums Sk{zi, . . . , z n ) = Y0£=i z \- I n their work M.Stone et 
al. conjecture that these polynomials correspond to the bosonic creation operators of edge 
modes. 

In this section the problem is reformulated in the language of second quantization and 
the second quantized counterpart of Sk are shown to describe edge excitations as free bosons 
in the limit n — > oo. Next the dispersion relation due to the interaction potential is dis- 
cussed also relying on a Hartree approximation. We shall, as in the work of Stone et al., 
assume no confining potential so that the Hamiltonian consists only of the e~- e~ interaction 
energy. The competition between the confining potential and the pair interaction in the EW 
spectrum is discussed elsewhere@. While bulk excitation involve the small q part Nx(q) 
we expect that edge excitations depend on the large T behaviour of this distribution, since 
electrons on the edge of the quantum dot belong to pairs with T ~ 2n. On the scale of the 
total angular momentum L = L a + M, a simple hydrodynamical argument^ shows that edge 
excitations involve changes AL = M ~ y/n (contrary to quasi particle excitations for which 
AL — M ~ n). In the thermodynamic limit typical values of M and k are of the order ^Jn. 
Note also that the size of the Hilbert space depends only on M as long asM<n. 

Sk is, in second quantized form, the single particle ladder operator 

&k = V ~T~^~ C m+k C m (18) 

m=0 V m - 

The conjugate operators are easily defined. In general these do not satisfy commutation 
relations typical of bosonic creation and destruction operators. However M. Stoned showed, 
using polynomial algebra, that, in the limit n — > oo and for v ~ 1, they do form a bosonic set 
of creation operators in the sense that the overlap between states with different occupation 
numbers (i.e. with different combinations of S£) vanishes. The same can be shown to hold 
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in the framework of second quantization for the operators S k . If \k) = A k (n)S k \i/jo(n)) and 
\k,j) = Akj(n)Sfr Sj~\ijjo(n)) are normalized states, a straightforward calculation shows that 
Ak(n) oc rT k l 2 and (k,j\k+j) oc n -1 . The general m bosons state, in the sector L = L Q + M, 
is 



\k h ...,k m ) = A kl> ^ km (n)S^, SjJM 71 )) 



(19) 



where Y^ih = M and A kl) .„ jkm (n) oc n~ M l 2 . The overlap between states with different 
occupation numbers vanishes at least as n" 1 . The state \ij)o{ n )) is the vacuum state for the 
bosons created by S k since it is annihilated by S k for any k. 



If these were free bosons we would have 



H,Si 



fl k S£. This commutator, if the 



Hamiltonian H has the form fl3T), can be written as: 



H, S k 



T=l q 

+ V T (q) [ft(q),S+\f T {q)} 



+ 



If C k i(T, q) is defined such that 



f T+k (q),S+}=J2C k , l (T,q)f T (q-l). 



(20) 



then0 [/t (9)5 S k = Ck,i(T, qjfr+kiQ + 0- The coefficients C k j(T, q) are easily evalu- 



ated with the help of equation (|7|): 



C k ,(T, g) = 4 ]T br{u, q)J { -^±^b T+k (u + k,q + l) = 2(-l)< (* 



u=0 



{T + k-q-l)\{q + l)\ 



l)\ 2 k (T-q)\q\ 



(21) 



clearly I < k otherwise C kl (T,q) = 0. Since only even values of / occur, we will neglect 
factors (— I) 1 and implicitly assume, from now on, that J2i runs only on even values. The 
final outcome reads: 



St 



J2 C k j(T, q) [V T+k (q + I) - V T (q)] f+ +k (q + l)f T (q) 
T,q,l 



(22) 
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that is: S£ promotes a pair in the state (T, q) to states of larger TAM and RAM (T+k, q+l)- 
Ck t i(T,q) describes how k additional units of angular momentum are distributed among 
the relative and the center of mass motion of the pair f^ +k (q + I). The energy, for any 
such transition, changes by an amount Vr+k(q + — Vr{q)- Since C k=ll (T,q) oc <^o> if 
Vr(q) = e g does not depend on T, we find H, S£ = 0. This is a known resulti!10: 
Si create zero energy excitations because these concern translations of the center of mass. 
Another consequence of equation fl22"|) is that in the case of the HC potential the commutator 
turns out to be proportional to /r+A^l + 0/t(1) so that it annihilates the ground state for 
v < 1/3 thus yielding Q k = for all fc. Note finally that the EW spectrum is unaffected by 
a constant shift of the pair potential Vr(q) — > Vo + Vr(?)- 

The operator form of equation (p^) strongly resembles the action of a ladder operator in 
the relative and center of mass coordinates. Indeed the same procedure of section I leads, 
for the first quantized operator Ak,i(zi,Z2) = Z k ~ l £}, to the following expression: 



(T + k . q . mq+]l )Mq) 



(T-qW 



This is easily expressed in S k operators by expanding Z k ~ l £} in the individual particle 
coordinates Z\ and z 2 - The second quantization procedure applied to the resulting expression 
reads 



k 

k 



9=0 



-^(l + x) k - l (l -xY 
9 ] - 



J x=0 

2 1+ HoS+ (23) 



Equation (|22|) can be expressed in terms of operators A kl provided that Vr+k(q + I) — 
Vr(q) is independent of T and q. Actually this is only fulfilled by Vr{q) = — 7<? which 
corresponds, in real space (see eq. (|TT|) ) , to an harmonic potential v(r) = 27(1 — r 2 ) in 
the inter-particle distance r. This potential has been studiedl^E in first quantization and 
is a rather unphysical pair interaction being unbounded from below as r — > 00. It has 
been pointed out recently^ that this harmonic interaction leads to the disappearance of the 
FQHE. In the present context the stability of the FQH ground state, that depends on the 
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competition between the pair interaction and the confining potential, is assumed. Given the 
harmonic potential e q = — r yq, equation ( p2"D is readily translated into 



_> 



"27 Eft <4E* 



\2-^b k (l,g)jS+_ g S+ 

here the sum on I runs only on even values. The second term of (p3|) does not contribute 



(this is true whenever V T (q) = e q ). Since 2Zy / (jJ2- fc - 1 = fcfc fe (£, 0) + Vkb k (l, 1), using (§) the 
sum in braces gives k(5 g fl + 5 9j k) — Vk(3 g ,i + 5 g ^-i) and finally 

H, S£] = -2 7 A; (nSjf - (24) 

Note that this is an exact result (Sq = n). When evaluating the Hamiltonian for L = L Q +M 
using the basis set ( ]T9|) the first term will contribute to diagonal elements while the second 
to off-diagonal ones. The larger off diagonal element is the one between states | . . . , k, . . .) 
and | . . . , k — 1, 1, . . .). The magnitude of the latter will be of order n' 1 with respect to 



diagonal elements because of the explicit factor of n in equation (0). Other off diagonal 
elements, being proportional to the overlap of states with different occupation numbers, will 
be at least of order n~ 2 with respect to diagonal elements. 

In the language of standard perturbation theory the first term of equation (|24l) may be 
regarded as coming from the unperturbed hamiltonian H , the second from a perturbation 
XV. The validity of perturbation theory depends on the ratio between A and the separation 
AE between unperturbed eigen-energies (diagonal elements). Equation (|24f) immediately 
yields A/ AE ~ k/n. First order perturbation theory is then exact as n — > oo even for 
physical wave vectors for which k oc yfn: S k create free edge excitations with Q k = 2wyk. 
The generalization to more physical pair potentials insists on this same argument. Before 
turning to the general case it is useful to remark that for k = 1 the correct result £lk=i — is 
recovered and that the resulting EW dispersion is linear. The spectrum is then degenerate 
since £lk+m = &>k + ^m- If 7 > the boson energy is an increasing function of k. In general 
if e q > is monotonically decreasing Qk is monotonically increasing in k. 
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An approximation scheme has to be introduced at this point to deal with more physical 
interaction potentials. The approximation essentially consists in taking out of the sum on T 
and q an effective value Jk,i(n) of Vr(q) — Vr+k(q + 0- The T dependence of Vr(q), together 
with the second term of equation (|23|), will be dropped from now on. In practice the effective 
value of jk,i( n ) ma Y be evaluated by evaluating the / th term of equation ( |2"2] ) on the ground 
state: 



Z k,l i ,, 



(25) 



where Z k j is such that J2 q Pk,i(q, n ) = 1- The evaluation of Pk,i{q, n) is a tedious task of 
algebra that is omitted here since it leads to a complex and lengthy expression from which 
it is hard to extract the interesting properties^. 

With this approximation equation ( p2"D can be expanded in A kl operators and finally in 

operators 

H,S+] = -2£ (^jjk^Ati = ^E^, 9 H^ + -A + ( 26 ) 



where f4 j9 (n) = 4n 



'k\ lk,i{ n ) 



f)9 

9 ] - 



(27) 

J x=0 



The same considerations following equation (^4j) show that only the g — 0, k terms are dom- 
inant as n — > oo. Moreover in this limit also the approximation e q — e q+ i ~ 7fc,z(n) becomes 
exact for smooth monotonic potentials. Here smooth means that e q — e q+ i introduces a neg- 
ligible dependence on q compared to the dependence of Ck,i{T,q) on the same variable. If 
the scale of q values is proportional to n this condition is satisfied by very general potentials. 
The latter condition is verified if Pkj(q,n) depends on q only through the ratio q/n: 

Pk,i(q,n) = -^k,M/n) (28) 
n 

The reason for this is that the dependence on n comes into Pk,i{q, n) through the pair 
distribution functions Nx{q) which has this property. Also one can argue that the typical 
value of q for a pair of electrons on the edge of the sample is of order n since the inter-pair 

13 



distance, y/q, should be of the order of the disk radius R(n) ~ \/n. Equation ( [28]) has been 
verified numerically. Figure |] shows a very good collapse of nPk,i(xn,n) for n = 30,40 and 
k = 8, I = 2,4,6,8. 

The condition on monotonicity is relevant because otherwise e q+ i — e q would change sign 
for some q. The operator structure of H, S k could in this case be different from that 
of S£_ g S^. Another source of troubles if e q is not monotonic is that the sign of ■jk,i( n ) 
may change for different I and this may eventually cause cancellation in the g — 0, k term. 
Instead, for monotonic potentials e q+ i — e q , and thus r )k,i{ n )i has always the same sign, so 
that the g = and g = k terms in fl27|) are always at least of the same order of magnitude of 
the other ones. For example table | lists the overlap between the state (k\ and the ground 
state in the sector L = L Q + k and n = 30, for the Coulomb, the HC potential and for the 
potential e q = q/(q 2 + 16). The latter is not monotonic and has a maximum at q = 4. We 
see that in the former cases the overlap is very close to unity, while for e x this is not true for 
k > 5. At k = 6 the state (k\ is between the ground state and the first excited state while 
for k — 7 it is very close to the second excited state. Another evidence of the correctness of 
the approximation is that the overlap depends weakly on the potential. 

We can safely conclude that S£ are the creation operators of the edge modes 

H, S£ = —flkfi(n)S^ as n — > oo 

where Qk,o — &>k,k is given by (Wn)- In the case k 3> 1, a rough estimate of Qk,o{n) is 
obtained observing that in (|27|) for g = 0, k a binomial average is performed of 7fe,z(n) and 
then Q kfi ~ 2nj k>k / 2 {n). 

Let us analyze in more detail the EW spectrum resulting from a general potential. Of 
particular interest is the behaviour of Qk,o{n) for a fixed k as n —>■ oo. This is related to the 
behaviour of the dispersion relation of EW in the long wavelength limit. In fact the physical 
wave vector k is proportional! to k/^/n. The hydrodynamical picture! assumes a dispersion 
relation linear in k for k —>■ 0. The velocity c of EW is given by 

c = lim^= lim v^oH (2Q) 
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As long as the potential e q is monotonic, 7fc,z( n ) will always contain a linear term in I so 
that the linear term in k of Qk,o will always be present. The evaluation of the explicit k 
dependence of flk,o(n) is complicated by the dependence introduced by Pk,i(q,n) which is 
difficult to analyse. The scaling property of Pk,i(q, n), equation fl28|), allows however to draw 
conclusion on the dependence on n of Q^^n). These conclusions will be further confirmed by 
an Hartree approximation for v = 1 that will also give some indication on the k dependence 
of the dispersion relation. 

The method is based on the separation of the dependence on I and q of e q — t q+ i. If e q — 
e q+ i = a(l)(3(q) then j k ,i{ n ) = a (l)Vk,i( n ) where r)k,i(n) = J2 q Pk,i(q,n)(3(q). Tne dependence 
on n of rjkj(n) can be extracted using equation (|28|) , i.e. 



and P(q) = e~ q ' q ° that in equation ( p0|) means that the behaviour of r]k,i(n) as n — > oo is 
related to the behaviour of fik,i{%) as x — > as could be expected since only pairs with small 
RAM contribute to the energy. Also note that for I ~ k/2 oc y/n all the dependence on k 
and n is in f]k,i{n) since a(k/2) ~ 1. However the analysis of ^k,i{x) as x — > is a very 
difficult task. We will circumvent this difficulty introducing the Hartree approximation in the 
following. A point worth of mention here is that the behaviour in n of flk^n) for short range 
potentials depends on the properties of the ground state, i.e. on No assumption 

on the form of fj>k,i{x) is instead necessary to extract the n dependence of Qk,o{ n ) i n the 
case of a long range potential of the form e q = q~ a (here a = 1/2 would correspond to the 
Coulomb potential). The separation of the / dependence from the q dependence is possible 
in this case using power expansion. The term l m will have a coefficient (3 m (q) oc g~ a ~ m 
that in equation ( p0|) will give a term of the order of n~ a ~ m . In the thermodynamic limit 
the m = 1 term is dominant and all the others can be neglected even when / cx <Jn. 
The situation is then very similar to that of the harmonic potential since e q — t q+ i can 
be replaced by lYj q Pk,i{q 1 n)aq^ a ^ 1 . The resulting spectrum Qk,o will be linear in k in a 




(30) 



In the case of a finite range potential of the form e q = e q l q ° we may take a(l) = 1 — e 
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first approximation. However, due to the additional dependence on I and k introduced by 
Pk t i(q,n), the EW frequency could contain also higher powers of k. The coefficient of the 
linear term in Qk,o{n) will be of the order of n~ a that in equation fl29|) yields a sound velocity 
c ~ n l / 2 - a that vanishes for a > 1/2. In real space this means that for potential vanishing 
faster than 1/r, that is faster than the Coulomb potential, the contribution of the interaction 
energy to the EW velocity is zero. 

Let us now turn to the Hartree approximation to the EW frequency. The derivation of 
this approximation is presented in the appendix. Here we briefly comment on its nature 
before turning to the discussion of the results. The final formula can also be obtained from 
equation ( p2|) by replacing the pair of operators + 0/t(?) with their commutator 

that has a form similar to S£ (see (|15|)) of a density excitation. The justification for this 
is that the neglected term, f^Q) fr+ki^ + 0> tries to create pairs with TAM T + k on the 
filled Landau level and this is of course not possible if T + k < 2(n — 1). Since the terms 
with T > 2(n — 1) of equation ( |22| ) vanish on the v = 1 state (since fj'{Q)\' l l J o{ n )) = 0) the 
neglected exchange term acts only on the " Fermi surface" . Apart from this, the derivation 
also makes use of equation and then the condition q -C n is assumed. The final result 
reads 

oo / u \ j 

= E b £^ } (g/n) (31) 

j = l \ n J q 

The first two terms j = 1 and 2 terms were worked out explicitly and their expansion in 
powers of y 1 q/n were found to contain only odd powers. The dependence on q through the 
ratio q/n stresses once again that edge wave excitations involve pairs of RAM of order n. 
Note that, for k oc y/n, the linear term in k dominates on all the others and this suggests 
that the asymptotic spectrum is linear. This means that the harmonic approximation for 
the Hartree potential is a good one0. 

Since T^(y) oc + 0(y 3 / 2 ), it is straightforward that for a finite range potential 
^k( n ) ~ kn~ 3 / 2 . This result coincides with the power law behaviour found by Stone et.al.i 
numerically. If combined with equation (j30|) , it suggests that Hk,i{x) ~ x 3 ^ 2 for x 1. Note 
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that for a short range interaction also the possibility of a finite dispersion relation with a 
long wavelength behaviour of the type (k/y/n) 3 is ruled out by equation ([31]). The above 
mentioned results for a long range potential are nicely recovered together with the numerical 
result! fifc(n) ~ krT 1 ! 2 for the Coulomb potential. Moreover equation ( pTj ) provides a further 
indication for the asymptotic linearity of the dispersion relation for long range potentials. 
Finally the coefficient of the k 2 term, that removes the degeneration of the linear term, were 
found to be positive so that, as conjectured by Stone et al.i, the ground state in the sector 
L = L a + M is actually S^\i(}o(n)) . 

In summary it has been shown that for monotonic potentials edge wave excitations near 
v — 1 are free bosons as n —>■ oo. The creation operators of edge excitations are S^. This 
result is in some sense a consequence of the completeness of the basis set created by S£. 
Every operator that raises the total angular momentum of the system can be expressed as 
a sum of combinations of S£. This has been explicitly done for the commutator [H,S£] 
which turn out to have a dominant contribution, as n — > oo, proportional to S£. 

The behaviour of the EW spectrum has been discussed for general classes of potentials. 
In particular the long wavelength limit, that is related to the n — > oo limit for k ~ \/n, has 
been explored in a general way using the fact that pairs with a RAM of order n are involved 
in edge excitations. This statement, that refers to the difference in the RAM distribution 
of the excited state with respect to that of the ground state, is also displayed in the results 
of the Hartree approximation, where again the dependence on the RAM q comes through 
the ratio q/n. The fact that macroscopic quantum numbers are involved in edge excitations 
supports the validity of the classical hydrodynamical pictured. The special role played by 
the harmonic e q = —'jq interaction and the result, from the Hartree approximation, that 
the asymptotic dispersion is linear in k, is also reminiscent of a classical elastic response. In 
some sense this potential comes out naturally in dealing with edge excitations for any pair 
interaction e q ; the effective 7(71) being some average of d q e q . The reason why this comes 
out is essentially the same for which the single particle confining potential, coming from 
Coulomb interaction with nuclei, is usually modeled by a harmonic one@00'lil'@. 
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A final consideration concerns the extension of these results to v = 1/m. Note that 
equation ([24]) for the potential e q = —jq is an exact result for any v. This formula, as 
equation (|26|), is a statement about operators. 

There are two basic conditions on the ground state \ipo(n)) that have been used: 

i) S£ provide an orthonormal basis of the Hilbert space for L = L + M. The overlap 

between different states (ji, . . . , j m \k\, . . . , k p ) vanishes as n _1 . 

ii) annihilate the ground state for any k: Sk\4>o(n)) = 0. 

If these two conditions are satisfied also by the Laughlin state then the whole description of 
EW presented for v = 1 can be applied at v = 1/m. Actually if condition i) holds, the second 
one follows by a simple argument. This is because if it were not so also would provide 
low lying excitations so that the ground state energy would not have a cusp at v = 1/m. 
Evidences for the validity of condition i) for v = 1/m come from the parton construction 



for the HC interaction. 



of the Laughlin states0i0 and also from the result 
These strongly suggest that the same picture outlined for v = 1 holds for v = 1/m. This 
supports the idea@ that edge states at v = 1/m are in a 1 to 1 correspondence to edge 
states at v — 1. Another consequence of this is that, while the asymptotic behaviour of the 
EW spectrum is expected to change for finite range potentials, as it does change for the HC 
potential, going from v = 1 to v = 1/m (here Qk = 0), the same behaviour is expected 
for long range potentials. This is because the power law behaviour of Qk,o{ n ) on n depends 
explicitly on the structure of the ground state only in the former case. The conclusion that 
the Coulomb interaction gives contribution to the EW velocity while potentials vanishing 
faster as r — > oo do not, can then be extended to v = 1/m. 
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APPENDIX: THE HARTREE APPROXIMATION 



The pair interaction hamiltonian can be approximated by a single particle potential for 
v = 1 using fr(q)f T (q) = fMfr (?) - [/r(?)> fr (<?)]• If only the part of f}(q)f T (q) acting 
on the lowest n — 1 orbitals is considered the first term can be neglected. The hamiltonian 
then becomes V" = J2 u <n e u( n )ci c u where the energy of the orbital u is given by 

e u (n) = ^e q 2b 2 T (u,q) 

the sums on T and q run on all the values consistent with u < n and T — u < n for the above 
mentioned restriction. This single particle potential yields indeed the exact value ([16]) of 
the v = 1 interaction energy. The commutator with S£ is easily carried out with the result 



-J2[e m (n) - e m+k (n)} 



(m + k)\ 
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This can be approximated by —fl k (n)S£ by taking out of the sum the value of e m (n) — 
e m+ k (n) for m = n — 1 on the " Fermi surface" . In other words fl" (n) is the gain in energy 
if one electron on the outer orbital is raised by k units of angular momentum. The result 
can be cast in the form Q k (n) = J2 q e q^k(q, n) where 

2(n-l) 

r k (q,n) = 2 [ b 'r( n - 1 ,q)- b T+k( n - 1 + k ,q) 

T=q 



This expression can be evaluated using the variable x — (n — 1 — T/2)/yT/2 and the 
asymptotic form of b T (u, q) equation (^). The sum on T can be changed into an integral in 
dx whose upper limit can be taken to by oo for q T oc n (this is also the condition for 
the validity of @). After some algebra Tk(q,n) takes the form 

J (y) \4> 2 q {\fny) - A (y, 9)4>l (VnC(y, g))] dy 



IT JO 
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where J(y)dy, A(y,g) and C(y,g) are the expressions of y2/(nT)dT, JT/(T + k) and 
{n — l + k/2 — T / 2) / n(T + k) /2 respectively in the variables y — xj y/n and g = k/n. Note 
that A(y, 0) = 1 and C(y, 0) = y. Next the integrand is expanded in powers of g and the 
terms containing derivatives of <fi 2 q (x) are integrated by parts. The final result is equation 
(j3T|) where 

The dependence on q/n of the above expression comes from the fact that <f) 2 (x) falls off 
rapidly for x > y/q. The cases j — 1 and 2 were worked out explicitly with the result 
Fi(y) = y~'ly 3 + j^V 5 ~ ^kv 7 + • • • and F 2 (y) = \y- §y 3 + j§^y 5 - §^y 7 + . . .. 
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The function P k j(q,n) has a complex expression that depends not only on N T (q) but on 
the full "pair density matrix" Gr(q,p) = {i J o(n)\ff(q)f T (p)\tlj (n)). 

The commutator of a single particle potential V± = J2 m -^mC„c m with S£ yields ckS£ only 
if E m = cm that is for a harmonic potential. For a general potential the commutator is 
very close to S£ under the same assumptions of monotonicity and smoothness discussed 
in the text for e q . An effective harmonic potential can be defined for v = 1 with c(n) ~ 

pCtTti(ll m E m \m=n—l ■ 

See A.H.MacDonald and M.D.Johnson, Phys. Rev. Lett. 70, 3107 (1993) and references 
therein. 
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FIGURES 

FIG. 1. Plot of the distribution of relative angular momentum, N(q) = J2t^t(q), for a 6 
fermion system at effective filling 1/3. Lt t is the total angular momentum, Tn ma x is the highest 
Landau orbital considered. Lines are drawn only for guiding eyes. The full dots (•) refers to the 
exact ground state, open dots (o) to the m = 3 Laughlin state while the squares (□) refers to the 
state built with three pair creation operators of highest relative angular momentum. 

FIG. 2. Plot of the distributions //fc,j(<z/n) = nPkj(q,n) for n = 30,40, k = 8 and all even 
values I = 2, 4, 6 and 8. 
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TABLES 



k (km (km (km 

2 0.9995 0.9995 0.9995 

3 0.9967 0.9967 0.9967 

4 0.9896 0.9905 0.9885 

5 0.9758 0.9795 0.9611 

6 0.9525 0.9624 0.6077 

7 0.9168 0.9379 0.0275 
7 - - 0.9718 



TABLE I. Overlap between the state (k\ and the ground states of the Coulomb (|Vo))> °f the 
hard core (IV'o)) an d of the potential e x q (\4>o)) defined in the text. In the last line is the first 
excited state of the latter potential for k = 7. 
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